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Abstract 

In the present paper, we introduce a multiple Ramanujan sum for arithmetic 
functions, which gives a multivariable extension of the generalized Ramanujan sum 
i -^h ' studied by D. R. Anderson and T. M. Apostol. We then find fundamental arithmetic 

properties of the multiple Ramanujan sum and study several types of Dirichlet series 
involving the multiple Ramanujan sum. As an application, we evaluate higher- 
dimensional determinants of higher-dimensional matrices, the entries of which are 
given by values of the multiple Ramanujan sum. 
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O 1 Introduction 

In 1918, S. Ramanujan [17] studied the sum 

c{k,n):= 22 e ( k , nl ) = 22 ^d^' ^'^ 

• i-H , ; ( mod fc ) d\ gcd(fc,n) 

^ _ gcd(f,fc) = l 

where, in the first sum, I runs over a reduced residue system modulo k with gcd(7, k) = 1, 
e(r,n) := exp(2iry/— ln/r) and \x is the Mobius function. The sum c(k,n) is called 
the Ramanujan sum (or the Ramanujan trigonometric sum) and is widely investigated in 
connection with, for example, even arithmetic functions [8, 9] and cyclotomic polynomials 
[14, 15]. See [13] for the arithmetic theory of the Ramanujan sum. Among several 
generalizations and variations of c(k,n), D. R. Anderson and T. M. Apostol [1] (see also 
[3]) considered the sum 

S Lg (k,n):= Yl fi^Wl (1.2) 

d\ gcd(fc,n) 

where / and g are arithmetic functions. Clearly, Sf g (k,n) extends the right-most expres- 
sion in formula (1.1) and hence gives a generalization of the Ramanujan sum. 

Motivated by the study of the above generalized Ramanujan sum, in the present paper, 
we examine the following type of multiple sum for arithmetic functions j\, . . . , f m +i] 

S fl ,..., fm+1 (ni,..., n m+ i) := V /i(^-)/ 2 (^) • • • fm(^-) fm+i(d m ) , 

, . ,, — , «1 "2 «m 

dj\ gcd(n 1 ,...,n J -_|_ 1 ) 
(j' = l,...,m) 
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where gcd(ni, . . . , Tij+i) is the greatest common divisor of ni, . . . , We call this a 

multiple Ramanujan sum for fi, . . . , f m+1 . Notice that the above expression gives the 
generalized Ramanujan sum (1.2) when m — 1 and, moreover, the Dirichlet convolution 
h * • ' • * fm+i of A, ... , fm+i in the 'diagonal case"; m = ■ ■ ■ = n m+1 . 

The present paper is organized as follows. In Section 2, we introduce a multiple Ra- 
manujan sum ^/^'' /l+l w ith positive integer parameters 71, . . . ,j m so that Sf lr __j m+1 = 
Sfo'"f r 1 and study its fundamental properties as a multivariable arithmetic function, 
such as the degeneracies and the multiplicativity (see [21] for the theory of multivariable 
arithmetic functions). Then, since S^'.'.'jl+l belongs to the class of even arithmetic func- 
tions (mod rii) as a function of rii-, ■ ■ ■ , n m+ i in the sense of Cohen [9], we calculate its 
finite Fourier expansion, which any even arithmetic function possesses. This expression is 
important with respect to Section 4. In Section 3, we study several types of Dirichlet series 
having coefficients that are given by S^ 1 ' "'f^\- We treat not only single- variable Dirichlet 
series, but also multivariable Dirichlet series. For instance, an analogue of the formula of 
J. M. Borwein and K. K. Choi [4], which contains the classical Ramanujan formula con- 
cerning the divisor function a a , is obtained. Section 4 is devoted to a higher-dimensional 
generalization of the so-called Smith determinant [19]. We evaluate higher-dimensional 
determinants (hyperdeterminants) of higher-dimensional matrices (hypermatrices) , the 
entries of which are given by values of ^/^'"/m+i' ^ n f & ct, we derive a hyperdeterminant 
formula for even multivariable arithmetic functions. This includes the results of P. J. 
McCarthy [12] and K. Bourque and S. Ligh [5] for the 2-dimensional case, that is, the 
usual determinant case, and partially of P. Haukkanen [10], for the higher-dimensional 
case. 

We use the following notations in the present paper. The set of natural numbers, 
the ring of rational integers, the field of real numbers and the field of complex numbers 
are denoted respectively as N, Z, K, and C. For m, . . . , n*. G N, gcd(ni, . . . , n k ) (resp. 
lcm(ni, . . . , nk)) represents for the greatest common divisor (resp. least common multiple) 
of rii, . . . , rifc. For \x\ * s the greatest integer not exceeding x. We denote the Mobius 

function as fi(n), the Euler totient function as <p(n), the power function as 5 x (n) := n x 
for x G C, and the identity element in the ring of arithmetic functions with respect to the 
Dirichlet convolution * as e(ri) := [^J = 1 if n = 1, and otherwise. Note that 5°*fi = e. 
Throughout the present paper, we consider a product (resp. a sum) over an empty set to 
always be equal to 1 (resp. 0). 



2 A multiple Ramanujan sum 
2.1 Preliminary: 7-convolutions 

Let A be the set of all complex-valued arithmetic functions / : M — > C We always 
understand that f(x) = if x N for / G A. The set of all multiplicative and completely 
multiplicative arithmetic functions are respectively denoted by M. and M, c '. Namely, 
/ G Ai (resp. / G Ai c ) means that f(mn) = f(m)f(n) for gcd(m, n) = 1 (resp. for all 
m,n G N). As usual, the product fg G A of /, g G A is defined by (fg)(n) := f(n)g(n). 
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Let 7 G N. We define the 7-convolution of /, g G A by 

*?/)(*) ^E/^M^)- 

In particular, *i = * stands for the usual Dirichlet convolution. Define the function 
a 7 G A by a 7 (ri) = 1 if there exists d G N such that n = d" 1 , and otherwise. Then, it is 
clear that g* 7 / = g™*f, where := a 7 g. Note that the product * 7 does not satisfy the 
commutativity or associativity properties unless 7 = 1. We therefore inductively define 
the 7 = (71, . . . , 7 m )-convolution for 71, . . . , j m G N of fx, . . . , f m+1 G A as 

(/m+l *7m • • • *7i A) W '■ = ((/m+1 * 7m ' ' ' * 7 2 A) * 71 /l) (n) 

= E A ("p ) /2 (^2 ) ■ • ■ /m ( ) fm+1 (<C) ■ 

Let L(s; /) := Y^m=i f( n ) n ~ s be the Dirichlet series attached to / G A and let <r(/) G 
RU {00} be the abscissa of absolute convergence of L(s; f). Define G A for 7 G N by 
/W(n) := /(nT). Then, it is clear that L( 7 s;/^) = L(s;/ W ) with tr(/W) < cr(/). For 
fx, ... , f m+ i G ^4, we here calculate the Dirichlet series attached to f m+ x * 7m • ■ • * 71 fx- 

Proposition 2.1. Suppose 70I71I72I • ■ ■ \j m with 70 = 1. TTien, we /icwe /or Re (s) > 
max{(r(/i),...,(j(/ m+ i)} 

m+l 

^(s; /m+i * 7m • ■ • * 71 A) = n ( 2A ) 

Proof. Assume d%?\ ■ ■ ■ |<i 7l |n. Then, since 71J72I ■ • • |7m, we can write d^ 3 = rfJ+Y^J+i f° r 
1 < j < m — 1 and n = &+<i 71 = • ■ • = kxk^ 1 • ■ ■ k^-^d 1 ™ with £+, k 2 , . . . , k m G N. We then 
write d m = k m+ i, and so 

00 

/m+i * 7m • ■ • * 71 /1) = £ h(h)f 2 (k?) . . . / m +i(fcr+i) (w • ■ ■ cur 

= L( s; / 1 )L( 7 i S ; / 2 (72> ) • ■ -L( 7mS ; 

This ends the proof. □ 

Remark 2.2. We cannot expect that L(s; f m +i * 7m ■ • • * 7l /1) is expressed as a product 
of a Dirichlet series such as (2.1) for general 7 = (71, . . . , 7 m ) G N m . 

In the next subsection, we will introduce a multiple Ramanujan sum for fx, ... , f m+ i G 
A, which gives the 7-convolution f m +i * 7m • ■ ■ * 7l /1 in the diagonal case (for more detail, 
see Proposition 2.3 (iv)). 

2.2 Definition of Sj and its basic properties 

Let 7 = ( 7l ,..., 7m ) G N m and / = (/ l5 . . . , f m+l ) G .4 m+1 . We define a multiple 
Ramanujan sum = S^ 1 ' '"f" 1 ^ of / with the parameter 7 by 

S}(n u . . . ,n m+1 ) := £ /i(^)/ 2 (|r) " ' ' /m(-^) 

dj 3 I gcd(ji 1 ,...,Ji J _ ) _ 1 ) 
(j' = l,...,m) 
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where the sum is taken over all m-tuples (di, . . . , d m ) G N m satisfying dJ J | gcd(ni, . . . , rij+i) 
for each 1 < j ' < m. Note that the summand vanishes unless d%?\ ■ ■ ■ \cCy\n\. We write 
Sf := Sf"' 1 ^ and understand that Sf(n) = f{n) when m — 0. These are the fundamental 
properties of S~J, which are obtained from the elementary properties of the gcd-function. 

Proposition 2.3. (i) For any 1 < j < m + 1, we have 

S](m, . . .,n m+1 ) = s-C-n.-./w-O K . . . ^ (2 . 2) 

/iv,/j-i,5 /j ,..., /m+l (-,n j+ i,...,n m+ i) 

(ii) For any 1 < j : < m + 1, we /iare 

5J(ni, . . . . . ,n m+ i) = fj(l) ■ ■ ■ f m+1 (l) ■ S^^ffjf' \n 1 , . . . ,n,-_i). 

In particular, we have Sj(l, n,2, ■ ■ ■ , — /i(l) ' " ' 

(hi) For any 1 < j < m, we have 



(71. -,7m) 

: ./j + lv./m 

5 /wJiT,/Jit,™!/^i ("i' • • • ' n i-i> S cd K-' ra i+2 ; ■ ■ ■ , fWi) (2-3) 



l(71v,7j-1.7j + lv.7m) 

fj-llfj+ 



and Sfcj^K . . . )M ) = Sfcj^W . . . ,n m ). 
(iv) Let rixlrij /or all 2 < j < m + 1 . Then, we have 

S](n u . . . ,n m+1 ) = (fm+i * 7m ■ ■ ■ * 7l /i)(ni). (2.4) 

□ 

Example 2.4. Let a = (a 1; . . . , a, n ) G C"\ Define a multiple divisor function o~2 with 
the parameter 7 by 

^(m,...,^)^ £ czr 1 ---^- (2.5) 

d^\-\dp-\m 

d/ I gcd(m,...,Tij + i) (j=l,...,m) 

This is a generalization of the usual divisor function a a (n) := ^2 d \ n d a ; a a {n) = a^\n,n). 
Then, it holds that &l = S] with ft = <J<«o+«i+-+<*j-i for 1 < j < rrz + 1 where 
a = 0. Similarly, for b = {pi, . . . ,6 m+1 ) G C m+1 , one can see that c~(ni, . . . ,n m+1 ) = 



1 1 S'J bi ^ m+1 (n 1 ,...,n m+1 ) where b := (o 2 - h, o 3 - 6 2 , • • • , b m+1 - b m ) G C m . Note 
that the sum Z2{a\, . . . ,a m ) := o~2(n, . . . ,n) is studied in [11] and called the multiple 
finite Riemann zeta function. From formula (2.1), if 70 1 71 1 72] • ■ ■ \ jm with 70 = 1, we have 
L(s; Z~ 1 (ai, a m )) = Y\™=x C(7i-i( s - a - • • ■ - a,-_i)) where ((s) = L(s; 5°) is the 
Riemann zeta function. 

Example 2.5. Let / G A. Then, the composition function / o gcd of / and the gcd- 
function can be expressed in terms of the multiple Ramanujan sum. Actually, from 
the degeneracy formula (2.3), we have (/ o gcd)(ni, . . . , n m+ i) = S/(ni, . . . ,n m+ i) with 
fi = 5°, f 2 = ■ ■ ■ = f m = £ and f m+1 = f * fi. 
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We next show the multiplicative property of the multiple Ramanujan sum SJ. Recall 
that an arithmetic function F(ni, . . . ,n k ) of /c-variables is called multiplicative if 

F(m 1 ni 1 . . . , m k n k ) = F{m ly . . . , m k ) ■ F(m, ...,n k ) 

for relatively prime fc-tuples (mi, . . . , m k ) e N k and (ni, . . . , n k ) G N k . Here, we say that 
(mi, . . . , m k ) and (ni, . . . , n k ) are relatively prime if gcd(m,j, rij) = 1 for all < i,j < k, and 
equivalently, gcd(ni ■ • -n k ,mi ■ ■ -m k ) = 1 (see [21]). In this case, we have the following 
Euler product expression 

F( ni ,...,n k ) = l[F(p^\...,p a ^), 

p 

where a p j = ord p rij for 1 < j < k. Note that this is a finite product since F(l, . . . , 1) = 1. 

Proposition 2.6. The function SJ is multiplicative if j\, . . . , f m +i £ »M. 

Proof. Assume (ni, . . . , n m+1 ) and (fci, . . . , /c m +i) are relatively prime. Then, for all 1 < 
j < m + 1, gcd(riifci, . . . , rijkj) = gcd(ni, . . . , rij) ■ gcd(fci, . . . , kf) and gcd(rii, . . . , nf) and 
gcd(A;i, . . . , kj) are relatively prime. Hence, the first assertion follows in the same manner 
as the that in the proof of Theorem 1 in [1]. □ 

2.3 Finite Fourier expansions 

An arithmetic function F(r; ni, . . . , n k ) of /c-variables rii, . . . , n k is called periodic (mod r) 
if F(r; ri\, . . . , n k ) = F(r; n[, . . . , n' k ) whenever rij = n'j (mod r) for all 1 < j < k (see [13] 
for the case of k — 1). It is well-known that F is periodic (mod r) if and only if it has an 
expression of the form of 

r 

F(r;ni, . . . ,n k ) = ^ a r (/i, . . . , / fc )e(r, n^i) • • • e(r, n fc / fc ) (2.6) 
ii >---j^fe=i 

and the coefficients a r (/i, . . . , 4) are uniquely determined by 

1 r 

a r (Zi, . . . ,l k ) = — ^2 F(r;n 1 ,...,n k )e(r,-n 1 l 1 )---e(r,-n k l k ). (2.7) 

ni,...,n k =l 

Moreover, F is called even (mod r) if F(r; m, . . . , n k ) = F(r; gcd(rii, r), . . . , gcd(nfc, r)). 
Note that F is periodic (mod r) if it is even (mod r). Then, as shown by E. Cohen [9], F 
is even (mod r) if and only if it has an expression of the form of 

F(r;ni, . . . ,n k ) = ^ a r (di, . . . , d k )c{d x , ni) ■ ■ ■ c(d k , n k ) (2.8) 
<fi,...A|j- 

with 

Or(di,...,d*) = ^ F(r;Si,...,8 k )c(^,^-) •■■c(£, (2.9) 

<Sl,...,tffc|r 

We call expressions (2.6) and (2.8) finite Fourier expansions and coefficients a r and a r 
finite Fourier coefficients of F. 
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By definition, the multiple Ramanujan sum Sj(rii, . . . ,n m+ i) is even (mod ni) as a 
function of m-variables n,2, ■ ■ ■ ,n m +i f° r an J n i £ N. Then, let us calculate the finite 
Fourier expansions of SJ. To do so, we add one parameter (or, a weight) to Sf. Let £ be 
an arithmetic function of m-variables. Set 

S}*(n u . . • , n m+1 ) = S^;y m if( nil . . . , n m+1 ) 

E • • • > O/ifiMl) • • • f m (-$£-)fn+i(d*r). 

7 . "i 1*2 Urn 

dj J I god(ni,...,n^ +1 ) 

y=i,...,m) 

Here, Sf (ni, . . . ,n m+ i) is again even (mod ni) as a function of n 2 , . . . ,n m+ i for any 
ni G N and Sj = Sj' lm where l m (di, . . . , d m ) = 1. We also write Sf := Sf"' . As 
discussed in Proposition 2.6, 5"^ can also be shown to be multiplicative if f±, . . . , f m +i 
and £ are multiplicative. 

For / = (/i, . . . , / m+1 ), we set */ = *Vm, • • • , ^/i)- Further, for n G N, we 

set 

m 

'61(4,..., <W := (n^(T^))-«(f---.T)- 

Then, we obtain the following theorem. Roughly speaking, the finite Fourier coefficients 
of SJ can again be written as a multiple Ramanujan sum. 

Theorem 2.7. Write the finite Fourier expansions of SJ as 



ni 



S^ni, . . . , n m+ i) = ^2 a }'ii • • • ' l m+i)e(ni, n 2 h) • ■ • e(ni, n m+ i/ m+ i) 

'2:— ,'ro+l=l 

= E a /iii( rf 2, • ■ • ,d m+ i)c(d 2 ,n 2 ) • • -e(d m+1 ,n m+ i). 

d 2 ,...,dm+i\ni 

Then, the finite Fourier coefficients aj '■ ni and ctf^ of Sf are respectively given as 

af ni (l 2 , l m+1 ) = n^slf (m, l m+1 , . . . , Z 2 ), (2.10) 
aj^, . . . , rf m+ i) = V^f 1 ("i. X^' " " " > T 1 )- (2- 11 ) 

Proof. We only verify formula (2.11) (formula (2.10) can be similarly obtained). From 
formula (2.9), ^'^(ok, . . . , d m+ i) is given as 

nf £ Sj'Hni, 52,..., 6 m+1 )c(f, ■ • ■ c(-^-, -^-) . 
t , c 2 (X2 o m+ i a m+ i 

02,-..,Om+l|«l 

Further, by changing the order of the summation, this expression is equivalent to 



X 



e 1 ,1 |5 2 |ni e™|5 m+ i|m 
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Here, we use the following identity. Let e\n and d\n. Then, we have 



- if d|e, 
e|5|n ^ otherwise. 

Actually, one can obtain this formula from the right-most expression (1.1) of the Ramanu- 
jan sum c(k, n) and formula 5° * \i = e. Then, applying formula (2.13), we see that (2.12) 
can be written as 



n7 m 



EE'" E £(e?,..,e-) 



X /l(77r)^ 2 (7^") ' ' ■ fm{-^r) fm+l{e^) 

02 Cm t>i Cm 

Changing variables e' fc = ni/e k /k for all 1 < k < rn, we have a 7k (ni/e' k ) = 1. Moreover, 
it holds that e' fc |^-j- for 1 < k < m and e[ \ ■ ■ ■ \e' m \ni, because e^ n | • • ■ {ej 1 \n\. Hence, we 
can rewrite the above expression as 



E E -E<m(?)-^(?)-f(?.-.^) 

In-, e' le' e'le' 1 

I™1 e m _ll e m e ll e 2 

m+1 e m-l'd^ C l 1 "^i" 



e m\ n l e m _xl e m e ll e 2 



x /i /2 ($)■•• /m (;rS W (5-) ^ ■•■ e' m . 

e l e m-l 



' e 



It is easy to see that (2.14) coincides with the right-hand side of formula (2.11). This 
completes the proof of theorem. □ 

Remark 2.8. Suppose that To I Ti I ' ' ' 1 7m with 7 = 1. Then, we have 5" t ~ ni = S where 

*£ ni : = and *J H = (5° f [ ^\, 5 1 f^ m ~ l] , . . . ,5 m f[ lo] ). In fact, since the condition 

above means that a lk (ni/e' k ) = o,y k (e' k+1 / e' k ) for all 1 < k < m — 1, the summand in (2.14) 
is written as 

. . . , ^(^(f*) • • • /IT 11 (^)fti(^)e[ ■■■e' m . 



e ™ e l e m-l e 



Hence, the claim follows. 



Example 2.9. Retaining the notation in Example 2.4, one can easily see that */ = 
( ( 5 ai+ -+ a '", (5 ai+ -+ a ™- 1+1 , . . . , 5 ai+m - 1 , <5 m ). Hence, the finite Fourier coefficient a f>ni : = 
a f,ni ™ °f ^ ne mu hiple divisor function a a := Oa ' " is given by 

/ 7 7 \ oiH ha m -m / 71 1 

a/,ni(a 2 , • • • ,fl m +i) = V ^i-*o(rii,- ,...,—), 

«m+l «2 

where 1 := (1, . . . , 1) e C m and *a := (a m , . . . , a x ) E C m 

Example 2.10. Retaining the notation in Example 2.5, we have Sf = /ogcd. Again, by 
the degeneracy formula of (2.3), the finite Fourier coefficient a^ ni of / o gcd is given by 

a fm (d 2 , ■ ■ .,d m+1 ) = n7 m S' / ^ i 5-(ni,gcd(^-, . . . , -p—))- 

Q>2 "m+1 



8 



Yoshinori Yamasaki 



3 Dirichlet series attached to Sj 

In this section, we examine both single-variable and multivariable Dirichlet series, the 
coefficients of which are given by the multiple Ramanujan sum SJ. 

3.1 Single-variable Dirichlet series 

We first examine a single- variable Dirichlet series. Recall the following well-known formula 
concerning the Ramanujan sum c(k,n) (see, e.g., [20]): 

£ c (fc jn )*- = ^W (Re( S )>l). (3.1) 
In this subsection, we give a generalization of this formula. For j = 1, 2, . . . , m + 1, let 

OO 

$J(s; hj) := ^ S f( n i, • • • > n m+1 )nj s , 

where fij := (ni, . . . , n 3 -_i, Tij+i, . . . , n m +i) e N m . The following proposition says that the 
series $^(s; hj) is written as the product of a Dirichlet series and a finite sum, which is 
again given by the multiple Ramanujan sum. 

Proposition 3.1. (i) For j = 1, we have 

$}(s; n x ) = L(s; h)S Fl (n 2 ) = L(s; faFfa) (Re (s) > a(/ 1 )), (3.2) 

where F l = F?$ M := 8" (tffcwl ( " , "3, • • • , ?Wi) * 71 5 s ) . 
(ii) For 2 < j < m + 1, we /iai>e 

$}( S ;n,) = C(«)^;; > >^(n ll ... I n i _ 1 ) (He (a) > 1), (3.3) 

Proof. Suppose that g?^ I g c d(^i, • • • , fifc+i) f° r all A; = 1, . . . , m. Then, rij is a multiple of 
lcm«, . ..,d*?) = dj 1 if j = 1 and lcm^r/, • • • , < m ) = d] 3 ~l if 2 < j < m + 1 (note 
that we only consider m-tuples (di, . . . , <i m ) such that d^ \ ■ ■ ■ IdJ 1 ). Hence, for j = 1 and 
Re (s) > $^(s; Wj) is expressed as 

E Em|t)m4) •••/„(%£■)/-»(«•) WO"* 

-y, ;_i U l U 2 Um 

d k K \ gcd(n 2 ,...,n fe + 1 ) (-1 

(fe=l,...,m) 

and, for 2 < j < m + 1 and Re (s) > 1, 

n d 11 d rj ~ 3 d Jj ~ 2 

Yl /i(^n)M^y ■ ■ ■ fj- 2 {Jy-i) fi^ijfJ^) 

d 7 *Mgcd(n 1 ,...,» fe+1 ) 1 2 J" 2 d^ 1 !^" 2 j ~ 1 

(fc=l,2,...,j-2) 

E EA(if)---/4^)/™ + i(e)(^i 1 <r- 

d^lgcdfdji^.n^! n fc + 1 ) 2-1 J ' 

(i=j,..,m) 
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Thus, it is easy to see that these coincide, respectively, with (3.2) and (3.3). This com- 
pletes the proof. □ 

Example 3.2. For small m, the series $^(s; hj) is explicitly given as follows. For m = 1, 

we have 

na) = L(s-J\) E M**)*^ 1 ', (3.4) 

d^i \n 2 

= C(«) E /i(^)/^ 7l ) rf ~ 7 " (3-5) 

<£ti |m 

and for m = 2, we have 



J7l 

$^ 3 ( S ;n 2 ,n 3 )=L( S ;/ 1 ) £ E M^)^?) V, 



^Ina^lgcd^.ns) 2 



71 



4Sl(-;-^3) = coo E E /i(^)/ 2 (^)/3(^ 2 )V ls , 

d{ 1 \n 1 dl 2 \gca{d[ 1 ,n z ) 1 2 



71 



72* 



d~l 1 \ gcd(ni,n2) 

Formulas (3.4) and (3.5) were obtained in [1] for the case of 71 = 1. 

Example 3.3. For < k < m + 1 and a = (a±, . . . , a m+1 ^ k ) E C m+1_fc , we define 
generalizations of the classical Ramanujan sum 0(72,1,712) as 

C m+l,fc( W l> • • • 5 n m+l) '■= 5Vi, . . . , / U,5 a i,...,<5 a »"+i-fe \ n l-> • • • ) n m+l)- 

It is clear that 0(711,712) = c^\(nx, 712)- Then, for 1 < k < m, using formulas (3.2) and 
L(s; /i) = 1/C(s), we can verify by induction on m and k that 



E C+u-K, • ■ • , n m+1 ) ni Sl ---n k Sk 

m,...,nk=l 



n^CC%) ^-^+-+^) (7a (d,n fe+ 2,...,n m+ i), (3.6) 

nJ=iC(«i+...+ ai )^ 



where a := (a 2 — 01,03 — a 2 , ■ ■ ■ , a m+ i_fc — a m _fc) G C m ~ fc (see Example 2.4). Setting 
m = 1, = 1 (in this case, 05 = 1) and ai = 1, we obtain formula (3.1). See [6] for a 
multivariable analogue of formula (3.1). 

Example 3.4. Retaining the notation in Example 2.5, we have Sf = f o gcd. Let 
gcd(nj) := gcd(ni, . . . , n 3 -_i, ■ ■ ■ , n m +i). Then, since it holds that gcd(ni, . . . , n m+ i) = 
gcd(rij, gcd(nj)), we have the following well-known formula 

00 

E(/°g cd )(^i ; - • -,n m+1 )nj s = ^ s o )f ^(s;gcd(nj)) = C(s) E U * ^)W~ S ■ 

nj=l d\gcd(fij) 
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3.2 Multivariable Dirichlet series 

For an arithmetic function F{n\, . . . , n^) of /c-variables, we denote L(s; F) by the multi- 
variable Dirichlet series attached to F, that is, 



" s l . . . n~ Sk 



L(s;F):= ^ F ( n i> ■ ■ ■ > n k) n i 

ni,...,n k =l 

for s = (s l7 ...,s k ) G E(F), where S(F) C C fc is the region of absolute convergence for 
L(s; F). The first goal of this subsection is to calculate L(s; S^) explicitly. 

Theorem 3.5. Let 'To I Ti I T2 1 • ■ ■ 1 7m with 70 G N. Define the function a 10 Sf of (m+ 1)- 
variables as (a 70 5'J)(ni, . . . ,n m+ i) := a J0 (rii)Sf(rii, . . . ,n m+ i). Then, we have 

m+l m+l 

L(s; a 10 S}) = J] C(^) ■ ]J L(s x + ••• + *,•; f^ ] ) (3.7) 
i=2 j=i 

in i/ie region 

' Re(sj) > 1 (2 < j < m+l), 

Re (si + • • • + Sj) > o-(fj) (l<j<m + l) 

In particular, setting 70 = 1, we have 

m+l m+l 

l ^ s d = n c(*i) • /1) n + • • • + fi jj - i] )- (s-s) 

3=2 3=2 

Proof. This is proven by induction on m. Suppose that m — 1. Then, it holds from 
formula (3.5) that 



s — (si, . . . , s m+ i) G C m+1 



L(( Sl , S2 );a 70 ^ 2 ) = E $^ 2 ( S2 ;n? )ni 



70S 
ni=l 



00 

7l s 2„-70Si 



c(* 2 )£ E /OM^K'^ 



Here, since 70I71, n 70 is expressed as n]° = a" yi P° with Z G N. Hence, this expression is 
written as 

c(- a ) EE/' (^)M rf71 K 7iS2 (^^)- si , 



d=l Z=l 



thus completing the proof for m — 1. Next, suppose that m— 1. Note that Sl(ni, . . . , n m+ i 
0(71) 

/l ' 5 / a ....,Wi( , ' ns '- ,n 
L(s; a 7o Sf) is equal to 



S 1 *- 71 ^ j (^1,^2) from formula (2.2). Then, from (3.7) for m = 1, 

/l ' S /2.-'-.Wl ( '' n3 '---' nm+l) 



n 3 ' ' ' n m+l 



n 3 ,...,n m+1 =l ■"• ■'a.-.Jm+l 

00 

= a^)L( Si; ff o] ) E + ^ ^fc&lc ■ >"*■■■> n m+1 )K 81 • • • 

rt3,...,n m+ i=l 

= C(-a)^(«i; /l 7ol )^((^i + S3, ... , Wi); an^to^l). 
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This completes the proof for m based on the assumption of induction. Hence, we obtain 
the desired formula. □ 

Remark that formula (2.1) is regarded as the "diagonally-summed version" of (3.8) 



since 



Example 3.6. Retaining the notation in Example 2.4, we have 

771+1 

£(«; o"a) = C(si) n C(si)C(«i + • ■ • + Sj - ai a,-_i). (3.9) 

fc=2 

Example 3.7. Retaining the notation in Example 3.3, we have 

T( a \ _ ns 1 c(si) ixi 1 ^ c(fi + ■ • • + s k+ i - at) 

1j{S, c m+lk ) - k 

1L=1 C(Sl + • • • + 8i) 

We can also obtain this formula from (3.6) and (3.9). In particular, we have L((si, s 2 ); c) = 

C(S2)C(S1 + S2-1)/C(S1)- 

Example 3.8. Retaining the notation in Example 2.5, we have Sf = f o gcd. In this 
case, it holds that 

r / /■ i\ C( S l) ' ' ' C( S m+l) r / 

L(s; / o gcd) = — ■ ■ -L( Sl + ■■■ + s m+1 ; /). 

C(«i H 1- -Wl) 

Next, we present another type of multivariable Dirichlet series involving Sj. Recall 
the well-known Ramanujan formula for the divisor function o~ a (see e.g., [20], also [7]): 

C(s)as-a)as-b)gs-a-b) 

L{s; o a a h ) = — 

C(2s - a - b) 

for Re (s) > max{l, Re (a) + 1, Re (b) + 1, Re (a + b) + 1}. This formula was generalized 
by J. M. Borwein and K. K. Choi [4] as follows. Let fi, f 2 ,gi,g 2 £ A4 C . Then, we have 

£ («;(W0ta *,»)) = Z^TKto^H < 310 > 

for Re(s) > max{<r(/i^i), a(hgi), cr(fig 2 ), <r{h92), cr(fif 2 g 1 g 2 )/2}. They gave formula 
(3.10) for the case in which 7 = 1, and it is easy to obtain the equation for general 7 G N 
in the same manner. Similarly, since Sfl\ 2 ( n i n ) = (f2 *7 fi)( n )> we can regard formula 
(3.10) as the diagonally-summed version of the following formula. 

Theorem 3.9. Let fx, f 2 , g±, g 2 G M. c ■ Then, we have 

L({s,s 2 )-S^ h S^ 92 ) 

C(a 2 )£(si; hgi)L(s x + s 2 ; (/ 2 ^)W)L( Sl + s 2 ; (fig 2 p)L( 8l + s 2 ; (f 2 g 2 p) 



L(2( 5l + s 2 );(/ 1 / 2 ^ 2 )W) 

m i/ie region 



(3.11) 



:si,s 2 ) G C 2 



Re(s 2 ) > 1, Re(si) > a(/^i), 

Re (si + s 2 ) > max{cr(/ 2 5fi), a{f x g 2 ), a(f 2 g 2 ), cr(f 1 f 2 g 1 g 2 )/2} 
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Proof. First, note that, for a multiplicative arithmetic function F of ^-variables, L(s; F) 
has the Euler product expression L(s; F) = Y[ p Y^h i -1 ^{P ll i • • • ,p lk )p~ Sltl ■ ■ -p~ Sklk . 
Let hi, h 2 G M. c . Then, it holds that 

where m 7 (Zi, l 2 ) ■— niin.{/i, ^2}] +1. Hence, from the Euler product expression, we have 

L^s^S^All) (3-12) 

00 

=n e ^/ 2 (p zi >^ 2 )^? 92 (^ i '^ 2 k si v s2 ^ 

P h,h=0 

00 , 



Jl,i2=0 



where a* := / 4 (p), ^ := 4 (p), 2* := p * ! for i = 1,2 and C 7 = C 7 (/i, f 2 , gi, g%) ■ = 
(aj — a 2 )(bj — b 2 ) / '(aifci) 7 for each prime p. Write I for the inner sum of the right-most 
hand side of (3.12). Moreover, divide I into two parts as Ylhh=o = ^2h>i 2 +2\^ 1 <z 2 = 
Yli°=o Y^h=h ^«T=o Yli°=h+i an( ^ denote I\ and I 2 by the former and latter sums, re- 
spectively. Writing l 2 as l 2 = 7/ + k in l\ and l\ as Zi = 7Z + k in 7 2 with / e Z> and 
< k < 7 — 1, respectively, we have 

7— 1 00 00 , 
',=EE E (a 1 6 1 x 1 )"x2«(l-(^)* 1 ')(l-(|)* 1 ») 



fc=0 i=0 «i=7i+A: 

1 1 — (ai6iXiX 2 ) 7 
1 — ai^iXi 1 — aibiXiX 2 



M 



and similarly 



x 2 1 - (ai6ixix 2 ) 7 

i 2 — M 1 

1 — x 2 1 — aib\Xi%2 

where M := E^oK^i^) 7 ' (l - ) 7(m) )(l - (^) 7{ ' +1) )- Hence, we have 

1 - (athx^y 
1 — Ji + i 2 — 7- ry- rM. 

(1 - x 2 )(l - ai6ixi) 

Now, M is straightforwardly calculated as 

C^l-ja^hhxlxiy) 

(1 - (ai6ixix 2 ) 7 ) (1 - (a 2 bix 1 x 2 p) (1 - (a!6 2 XiX 2 )T) (l - {a 2 b 2 x 1 x 2 )~f) ' 

where the desired formula immediately follows from formula (3.12). □ 

Remark 3.10. For m > 2, we cannot expect that the single- variable Dirichlet series 
L(s; (f m +i * lm • ' • * 7 i /i)(flW+i * 7m • • • *7i Pi)) has a product expression such as (3.10) 
and, similarly, the multivariable Dirichlet series . . . , s m +i); ^/^'' /l+l^v'-im+i) has 

a product expression such as (3.11). 
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4 A generalization of the Smith determinant 

In this section, we evaluate the hyperdeterminants in the sense of Cayley of the hyper- 
matrices, the entries of which are given by values of the multiple Ramanujan sum Sj^. 

4.1 Hyperdeterminant 

Recall the definition and some properties of the hyperdeterminants in the sense of Cayley. 
Let A = (A(ji, . . . , *fe))i<ii,...,i fe <n be a A;- dimensional matrix of order n. For a subset 
I C {1, . . . , k}, set Ej — 1 if j G /, and otherwise. Then, the hyperdeterminant detjA 
of A with the signature / is defined as 

- k n 

' ai,...,cr k e<Bn j=l v=l 

where & n is the symmetric group of degree n. We obtain the usual determinant of a 
matrix A when k = 2 and / = {1,2}; det{ 12 }^4 = det A. Note that detjA is identically 
zero if the number of elements in I is odd. For the theory of hyperdeterminants, see, e.g., 
[16]. We next present some lemmas that will be needed in the subsequent discussion. 

Lemma 4.1. Let A = (A(z'i, . . . , i k )) be a k- dimensional matrix of order n. 

(i) For ii G & n , we have 

det/(A(7r(ii), . . . , n(i k )))= det I (A(i 1 , . . .,i k )). (4.2) 

(ii) For 7r G & k , se ~k 7r_1 (-^) := { 7r_1 (*) I i £ I}- Then, we have 

det/(v4(z 7r( i), . . . , «7r(fc)))= det^-i(/)(A(ii, . . . ,i k )). (4.3) 

Proof. These are immediately obtained from the definition of the hyperdeterminant. □ 

Let A = (A(ii, . . . ,i k )) and B = (B(ii, . . . be k- dimensional and /-dimensional 
matrices, respectively, of order n. Then, the Cayley product AB of A and B is a (k + l — 2)- 
dimensional matrix of order n given by 

n 

(AB) (ii, . . . , ik+1-2) ■■= ^2 y4 ^ 1 ' • • • ' *fc+J-2)- (4.4) 

j'=i 

Lemma 4.2. Let A = (A(z'i, . . . , i k )) and B = (-B(z'i, . . . , z';)) 6e k-dimensional and l- 

dimensional matrices, respectively, of order n. Let K C {1, 2, . . . , k — 1} and L C 
{2, 3, . . . , /} wzY/i o<id cardinality. Set I := K U (L + (fc — 2)), where L + {k — 2) : = 
{/ + — 2 | / G L}. Then, we have 

detj(AB) = det Ku{k} A • det { i }U L5. (4.5) 



Proof. See [10, 18]. 



□ 
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4.2 Smith hyperdeterminants 

Let S = {x\, . . . , x n } be a set of distinct positive integers. The goal of this subsection is to 
evaluate the hyperdeterminant deti(Sj'^(xi 1 , . . . , Xi m+1 )) for a factor-closed set S. Here, 
5* is called factor-closed if 5* contains every divisor of x for any x £ S. To accomplish this, 
we slightly extend the definition of the even arithmetic function examined in Section 2. 

First, we define some notations. For d = (d±, . . . , c4) £ N k , we set \d\ := X^=i dj and 
a^d := ■ • - ,x dk ). Moreover, for r £ N, we set gcd(d,r) := (gcd(di,r), . . . ,gcd(4, r)) £ 
write d|r if dj|r for 1 < j < k and, in this case, set r/d := (r/d±, . . . ,r/dk) £ N fc . 
Let f(m,n) be an arithmetic function of 2-variables. For m = (mi, . . . , m&) £ N k and 
n = (m, • • • , n k ) £ N fe , we set /(n, m) := nj=i %)■ 

Let r = (ri, . . . , r m ) £ N m and A; = (fci, . . . , fc m ) £ N m . Let = (7ij,i, . . . , fy-,^) £ 
N fcj for j = 1, . . . , m. We call an arithmetic function F(r; n 1; . . . ,n m ) of |fc|-variables 
ni, . . . , n m even (mod r^ fc ^) if it holds that 

F(r; gcd(n l9 r x ), gcd(n m , r m )) = F(r; n 1 , . . . , n m ). 

Then, the following lemma is obtained as described in [9] (note that the case of m = 1 is 
nothing more than Theorem 22 in [9]). 

Lemma 4.3. A function F(r; n±, . . . , n m ) is even (mod r^) if and only if it possesses 
a representation of the form 

m 

r \T , Tli, ■ ■ ■ i Tim) ^2 a r (d u . . . ,d m ) Y[c(dj,rij), (4.6) 

(j=i,...,??i) 

where the sum is taken over all dj = (c^i, . . • , £ N fcj for j = 1, . . . ,m satisfying 
dj\rj. Moreover, the finite Fourier coefficients a r (d±, . . . , d m ) are uniquely determined by 

aAd 1 ,...,d m ) = - K —— £ F(r;S l ,...,S rn )l[c(^), (4.7) 

1 m Sj\rj j=l 3 3 

(j = l....,m) 

where the sum is also taken over all dj = (5j,i, . . . , o~j,k 3 ) £ N kj for j — 1, . . . , m satisfying 
8j\rj. □ 

Now, we obtain the following proposition. 

Proposition 4.4. Let S = {x\, . . . , x n } be a factor-closed set, and let F(r; ni, . . . , «,&) 
be an even function (mod r^). Define the (m+ \k\) -dimensional matrix B of order n as 

B{i\, . . . , 2 m +|fc|) = Fyx^, . . . , Xi m , Xi m+1 , . . . , %i m+kl , • • • , ^i m+fcl+ ... +fcm _ 1+1 , • • • , x i m+ \ k \) ■ 



Let I be a subset of {1,2, ... , m+ |fe|} with even cardinality such that {m+1, m+2 . . . , m+ 
C I, and set Sj = 1 if j £ I , and otherwise. Define the subset I of {1,2, ... , m} as 
j £ I if and only if Ej + kj is odd for 1 < j < m. Then, we have 

detjB = (xi ■ ■ -x n ) detj(a Xiit __ iXim ( xj^ . . . ,Xj x , . . . , Xj m , . .j , a^ m J ) ^ m<n - (4-8) 

k i km 

Here, the finite Fourier coefficient of F given by (4.7). 
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We then need the following lemma. 

Lemma 4.5. Let A = (A(ii, . . . , i k )) and C = ((7(21,22)) be k- dimensional and 2- 
dimensional matrices, respectively, of order n. Fori = 0, 1, . . . , k, define the k-dimensional 
matrices Aq of order n by the following recursion formula: 

A^\i 1 ,...,i k ) :=A(i 1 ,...,i k ) (1 = 0), 

A<§ (z l9 . . . , i k ) := (A^C) (z 2 , ...,i k ,h) (/ = 1, 2, . . . , k), 

where A^ is the Cayley product of A^ and C . Then, we have 

n I 

A i §(i l ,...,i k ) = M i i+i,---, i k,ji,---,ji)Y[C(j h ,i h ). (4.9) 

ji,."!J/=l h=l 
Proof. This is shown by induction on / from the definition of the Cayley product (4.4). □ 

Proof of Proposition 4-4- From formula (4.2), we can assume that X\ < x 2 < ■ ■ ■ < x n . 
Define the (m + |fc|)-dimensional matrix A and 2-dimensional matrices C of order n, 
respectively, as 



A(ii, . . . , i m +|*|) : — a x n ,...,x im (xi m+1 , ■ ■ ■ , x i m+ \ h \) Y\ n e 
C(ix,i 2 ) : = c(x h ,x i2 ), 



j = llj=l 



where = 1 if Xj\xi, and otherwise. Let p = (1, 2, • • • ,m + \k\) G © m +|fc| be the cyclic 
permutation of order (m + |fe|). Then, we have 

B(ii, . . . , i m +\k\) = ^c fc '^(V m (i)' • • • ' V m (m+|fc|))- (4-10) 
In fact, by equation (4.9), the right-hand side of (4.10) is equal to 

A(j (i m _|_i, . . . , i m -\-\k\ 1 i\ 1 ■ ■ ■ 1 i>m) 

n 

^ ] ,—,3Hm (. X dl,l ' ' ' ' ' "^l.fci ' ' ' ' ' X dm,l 



' • • • ' •"1m.,kr, 

dj,l^<l j>k .=l V ' V - 

0=l,...,m) Kl 



x n n ed 3,ij' i j c ( xd 3,ij ' Xi m+k 1 +-+k j _ 1 +i j ) 



3=1 i 3 =i 



m Kj 

a ooi 1 ,-,x im ( x dii ■ ■ ■ i^dm) \\ \\ c { x d jJj ,%i m+kl+ ... +kj _ 1+l .), 



x d j \ x i j j=llj=l 
(j=l,...,m) 



where dj := (dj t i, . . . , dj tkj ) for 1 < j < m. This clearly coincides with the left-hand side 
of formula (4.10) from the expression of (4.6) of F, thus completing the proof. Note that, 
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since S is factor-closed, every divisor of Xj can again be written as Xd for some d. Next, 
using formula (4.3), we have 

detj.B = det/(A[l fc ^(i p m( 1 ), . . . , i p m( m+ \k\))) 

= det/((A^ fc ' 1 C) (i p m(2), • • • , i p ™( m +\k\), V m (i))) 
= det/((y4^ fc ' 1) C , )(i p m+i( 1 ), . . . , i p m+i( m+ |fc|))) 
= det p - m -i (/) ((A[i fchl) C , )(z 1 , . . . ,i m +\k\))- 

Here, note that m + 1 G / implies m + |fc| = p~ m ~ 1 (m + 1) G p~ m " 1 (I). Therefore, using 
formula (4.5) with K = p~ m (I) \ {m + |fe|} and L = {2}, we see that the right-hand 
side of the above equation is equivalent to 

det( p -m-i( / )\{ m+ | fc |}) U {, m+ | fc |}(74[l fc ' x) ) ■ det{i ]2 }C 

= detC • det p -™-i(7)((A^ fc| ~ 2) C)(i 2 , • • • , W|fc|> *i)) 
= detC ■ det p - m -2 (/) ((A^ fc|_2) C)(zi, . . . ,i m+ \ k \))- 
Under the condition m + 1, . . . , m + \k\ G /, the above procedure yields 

det f B = (det {lj2} C) |fc| ■ det p - m -| fe |(/)(^c 0) ) = (det {li2} C) |fc| ■ det 7 A 

By Corollary 3 in [5], we have detC = det(c(xi 1 , Xj 2 )) = x\ • • ■ x n . Hence, it is sufficient 
to calculate det/A Since Ej = 1 for j = m + 1, . . . , m + |fe|, we have 

<T (Tl,...,(T m j=l j=llj=l 

n m kj 

u=i i=i Zj=i 

Here, the sum is taken over all cr = (<ti, . . . , cr m ) G (6 n ) m and cr,, = (c^i, . . . , Oj^.) G 
(6 n ) fcj ' for 1 < j < m. For cr = (cr l5 . . . , a k ) G (& n ) k and t> G {l,...,n}, we set 
cr{y) := (ai(v), . . . ,&k( v ))- Replacing the variables r^\. = aj^.aj 1 with 1 < j < m and 
1 < lj < kj and writing Tj = (r J; i, . . . , Tj^.), we have 



<T j=l Tl,...,T m j=l L=l 



71 



v=l j = l l j= l 

where Tjdj = (t^iiTj, . . . , Tj^-ffj) for 1 < j < m. Note that e T . i a .^ v ) t(T ^ v ) = 1 holds for all 
1 < v < n if and only if r^y . = 1 since x\ < x 2 < ■ ■ ■ < x n . Therefore, detjA is equal to 



_^ m n 

—\ e n s g n (^-) £j+fcj n KiW' • • • ' 



<T j = l V = l 



— det J (a !til ,...,Xi m ( -ph, ■ ■ ■ ; ^iy • • • ; ^i m ; • - l X im) ) Kj^ j m < n " 

This completes the proof of the proposition. □ 
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Theorem 4.6. Let S = {xi, . . . , x n } be a factor- closed set and let F(r; rii, . . . , n^) be an 
even function (mod r) with respect to the variables rii, . . . , 7i&. Let I = {2, 3, . . . , k + 1} 
if k is even, and {1, 2, . . . , k + 1} otherwise. Then, we have 

n 

det I (F(x il ;x i2 , . . ■,Xi h+1 )) 1 < hr .. iim+1 < n = Oi • ■■Xn) k Y[u Xv (x v , . . .,x v ), (4.11) 

v=l 

where a Xv is the finite Fourier coefficient given by (2.9). 

Proof. Setting m = 1 in Proposition 4.4, we immediately obtain this theorem. □ 
Remark 4.7. The case of k — 1 is obtained in Theorem 2 in [5]. 

Corollary 4.8. Let S = {x\, . . . , x n } be a factor-closed set. Let I = {2, 3, . . . , m + 1} 

i/m even, and {1,2, . . . , m + 1} otherwise. Let £ 6e arbitrary arithmetic function of 
m-variables. Then, we have 

det 7 (S , J' ? (x il , . . ■,\ + i)) 1 < ilrijM+l <„ 

n 

= (A(i) • ••.unrn^''<- • • ■ ; ^)/K 71 '- ,7ro)] (^)- (4.i2) 

i id-7 

Proof. From formula (2.11), we have ay^(x^, . . . ,x v ) = x- m Si Xv {x v , 1, . . . , 1). Further- 
more, using the identity a yi (x v ) ■ ■ -a ym (x v ) = a\ cm ^ lt ,„^ m )(x v ), we can see that the right- 
hand side above is equivalent to x~ m ^(x v , . . . ,x v )fi(l) ■ ■ ■ f m (l)fm+i' Y1 ''''' ( x v)- From 
formula (4.11), we thus complete the proof. □ 

Example 4.9. Retaining the notation in Example 2.5, we have Sf = fo gcd. Hence, by 
formula (4.12), we have for a factor-closed set S = {x±, . . . , x n } 

detj((f o gcd)(x n , . . ■,lWi))l< ill ..^ 1 <n = 11^ * 

v=l 

This is a part of the result of P. Haukkanen [10]. 
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